DYNAMICAL SYSTEMS ADMITTING 
NORMAL SHIFT AND WAVE EQUATIONS. 



R. A. Sharipov 

Abstract. High frequency limit for most of wave phenomena is known as quasi- 
classical limit or ray optics limit. Propagation of waves in this limit is described in 
terms of wave fronts and rays. Wave front is a surface of constant phase whose points 
are moving along rays. As it appears, their motion can be described by Hamilton 
equations being special case for Newton's equations. In simplest cases (e. g. light 
propagation in isotropic non-homogeneous refracting media) the dynamics of points 
preserves orthogonality of wave front and rays. This property was generalized in 
the theory of dynamical systems admitting normal shift of hypersurfaces. In present 
paper inverse problem is studied, i. e. the problem of reconstructing wave equation 
corresponding wave front dynamics for which is described by Newtonian dynamical 
system admitting normal shift of hypersurfaces in Riemannian manifold. 



1. High frequency asymptotics for wave equations. 

Let M be a Riemannian manifold of the dimension n. Denote by V covariant 
differentiation with respect to metric connection in M . Then the expression 

eLr-D"- =J2-Y1 (-) ,2:")-Vfc, -...-Vfe, (1.1) 

fcl = l kr.=l ^ * ^ 

is an elementary differential operator of r-th order. Here a'^^'"'^^(x^, . . . ,x") are 
components of some smooth symmetric tensor field a of the type (r, 0) in local 
coordinates x^, ... , x" and i = \/— 1. Operators (1.1) are united into polynomial 

m 

H{p,D)^Y.''r-D\ (1.2) 



This is TO-th order scalar differential operator in M. Here p is a point of M and 
D is a formal symbol for differentiation. Operator (1.2) can be applied either to 
scalar field or tensorial field in M. We shall apply it to scalar field u, but first we 
introduce large parameter A in H{j),D). Let's denote 

m 

H{p,X-^D)^^\-''ar-D''. (1.3) 
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Then consider scalar m-th order wave equation with large parameter A: 

H{p,X-'^D)u = 0. (1.4) 

In particular cases (1.4) can reduce to Hclmholtz equation, to Schrcdinger equation, 
and to wave equation (see [1]). The following ansatz for formal asymptotic solution 
of the equation (1.4) as A — > oo was suggested by P. Debye: 

oo 

u = yJ^.e'^^. (1.5) 

Here ipt^a) and S are some smooth scalar fields in M. Scalar field S is interpreted 
as a phase of wave. 

Now let's substitute (1.5) into (1.4). This is done in several steps. First of all 
let's apply differential operator V^i • . • V^^ to e*^'^. This yields 

Vfe,...Vfe,e^^^= (^^/?<'^;...,^.(zA)-«j -e^^^. (1.6) 

Leading coefficient in this expression can be written explicitly 

/?r.fc. = Vfe.5-----Vfc.5. (1.7) 

Other coefficients f3'j?^ ^.^ can also be calculated for each particular q. It might be 
difficult to write general formula for all of them, but we do not need this formula. 

In the second step we apply operator (1.3) to power series (1.5). Using formula 
(1.6), for H{p, \~^D)u we obtain the following expression: 



oo m r r n 



If we collect all terms with fixed value of a and 6 = q + a,we find that the number 
of such terms in the above expression is finite. Simplest way to prove this fact is to 
make the following transformations of sums: 

oomrr oo m r q oo m m q 

EEEi:- 1:1:1:1: -1:1:1:1:' 

a=0 r=0 s=0 q=s a=0 r=0 q=0 s=0 a=0 13=0 r=q s=0 

00 ra m q 00 9 m 9 — a 

EEEE-E E E E- 

a=0 q=0 r=q 5=0 ^=0 a=aQ^ r=0—a s=0 



Here ae^m = max(0, 9 — m). Now we introduce differential operator R^^: 

Cl...fe,. fli^ 

s\ (r — s)\ 



m u n n rn} ^ki.-.k^ /Oi^^ — s) 

5-^ / . • « • Pi. , , ... 



^^ = EE E-E " \,r.-7r"^ -^^--^^-- (i-^) 



r=u> s=0 fci=l kr=i 
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Let's rewrite left hand side of the equation (1.4) in terms of operator (1.8): 

oo ^ D 

e=0 a=ae,„ ^ ' 

This is third step in substituting formal power series (1.5) into the equation (1.4). 
As a result wave equation (1.4) breaks into infinite series of equations for (fi(a)- 

6 

J2 R0-a^(a)=O for 6' = 0, 1, ... , 00. (1.9) 

In order to study the equations (1.9) we calculate explicit form of operator (1.8) 
for w = 0. Using formula (1.7), we derive 

m n n 

Ro = Y^y-y a''^-'"-{x\... ,x")-Vk,S-...-Vk^S. (1.10) 

r=Ofei = l kr=l 

Note that Rq is a scalar operator, i. e. it is the operator of multiplication by the 
function Rq = Ro{x^, . . . , x"). Therefore first equation in the series (1.9) looks like 
Rq ■ "^(0) = 0- Leading term in power series (1.5) for u is nonzero. Therefore the 
equation Rq ■ y>(0) = reduces to the following one: 

Ro = 0. (1.11) 

Taking into account (1.11), we can rewrite (1.9) as 

9-1 

Rg-a'P(a)=0, where ^ = 1, 2, . . . , oo. (1-12) 

First equation in the series (1.12) is also very simple. It is a partial differential 
equation of the first order for the function (^(q, : 

Rmo, = 0. (1.13) 

Other equations in the series (1.12) can be written in recursive form: 

e-2 

Ri'P(e-i) = - Ro-a<P(a) for 6* = 2, 3, . . . , oo. (1.14) 

«=«9,m 

Below we calculate explicit expression for operator Ri , but we do not consider the 
problem of solvability of equations (1.13) and (1-14). This should be done for each 
particular equation (1.4) taking into account initial-value and/or boundary- value 
problem stated for this equation. 

In order to find explicit expression for operator Ri we should first find explicit 
expression for coefficient /?^^' in (1.6). This is done by direct calculations. Ap- 
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plying operator Vfe^ . . . Vfe^ to e*'^'^ and collecting all terms with {i Xy ^, we get 

f3l\^ = for r = 1, (1.15) 

C.fe.= E Vfe,V,,.5.n V,,5 for r>2. (1.16) 

Now let's substitute the expressions (1.15) and (1.16) into the formula (1.8) and 
take into account symmetry of coefficients a''! This yields 

m n n r 

^1 = E E - E • """■■"'^ • n ^fc.^ • ^'^i + 

r=l fei = l kr. =1 p=2 

r{r — 1) 



m n n 

EE-E^^-v..v.,5.nv.,5. 

r=2A:i=l fer =1 P=3 



Further we assume that the equations (1.13) and (1.14) with operator (1.17) in left 
hand side are solvable in some sense, and we assume that that (^(qj ^ 0. Under 
these assumptions we concentrate our efforts in studying the equation (1.11). 



2. Hamilton- J ACOBi equation. 

Note that right hand side of the equality (1.10) is similar to that of (1.1). This is 
not casual circumstance. In order to reveal the nature of this analogy let's denote 
VS by p. Then p is a covector field with components 

Pk = ^kS, where k = 1, . . . , n. (2.1) 

Using notations (2.1), we can rewrite (1.10) as Rq — H{p,p), where 

m n n 

^f(p,p) = EE-E "'"••''^(^''•• - ^^'')-Pkr-----Pkr- (2.2) 

r=0/si=l kr=l 

Right hand side of (2.2) is a polynomial in components of covector p. If pi, ... , p„ 
are treated as independent variables representing components of some arbitrary 
covector p at the point p G M, then polynomial H{p, p) is called the symbol of 
differential operator (1.2). Once H(p,p) is given, the operator (1.2) itself can be 
obtained by substituting p^ = —iVk into its symbol. This trick is well-known in 
quantum mechanics, where momentum vector p is replaced by differential operator. 
The analogy with quantum mechanics becomes more transparent (see [2]) if we take 

Pk = -j^k (2.3) 

with A = l//i treated as large parameter. If we substitute (2.3) into the polynomial 
(2.2), then we obtain the operator (1.3) used in the equation (1.4). 

Now let's return back to the equation (1.11). It is first equation in the series 
(1.9). Due to (1.10) and (2.2) we can write it as follows: 



H{p,VS)=0. 



(2.4) 



DYNAMICAL SYSTEMS ADMITTING NORMAL SHIFT AND 



5 



The equation (2.4) is first order nonlinear partial differential equation with respect 
to scalar field S. It is often called Hamilton-Jacobi equation. Traditional way 
of solving such equations is based on the concept of characteristic lines. For the 
equation (2.4) these are parametric curves determined by Hamilton equations: 

In paper [3] the following form of Hamilton equations (2.5) was suggested: 

x' = V'H, VtPi = -ViH. (2.6) 

It was especially designed for the case of Hamiltonian dynamical systems in Rie- 
mannian manifolds. Here Vt is standard notation for covariant differentiation of 
vector-valued and tensor-valued functions with respect to parameter t along the 
curve. In our particular case covariant derivative WtVi is written as 

n n 

fe=i j=i 

Covariant derivatives Vj and V" in (2.6) are less standard objects. The matter is 
that, writing Hamilton equations (2.6), we treat H as extended scalar field in 
the sense of the following definition. 

Definition 2.1. Extended tensor field X of type (r, s) in M is a tensor- valued 
function that maps each point q G G C T*M to a tensor of the space T^{p,M), 
where p = 7r(g). Subset G of T*M is a domain of extended tensor field X. If 
G = T*M, then X is called global extended tensor field in M. 

In local coordinates extended tensor fields are represented by their components: 

^t't = . . . . . . ,p„). (2.8) 

In contrast to traditional tensor fields, components of extended tensor field are 
functions of double set of arguments. First n arguments of in (2.8) are 

local coordinates of a point p € M. Others are components of covector p € T*{M). 
Both p and p form a point q = {p, p) of cotangent bundle T*M. Further we shall 
not come deep into the theory of extended tensor fields, referring reader to paper 
[3] and thesis [4] for more details. But we give explicit formulas for covariant 

■ y X , Pi, . . . , Pn- 



dPq 



v«x;^\-;.; = (2.9) 



a=l 6=1 

r n s n 

E\ ' pifc -vil--- ak... ir _\^\ ' pfcfc vil 
Z^'^qak^ji is qjk h---bk- 

/s=lafc=l fc=l 6(5 = 1 



(2.10) 
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Covariant differentiation V defined by formula (2.9) is called momentum gradi- 
ent. Second covariant differentiation V defined by formula (2.10) is called spatial 
gradient. Note that for traditional tensor fields formula (2.10) reduces to standard 
formula for covariant derivatives. Therefore use of symbol V in section 1 and in 
formula (2.4) does not contradict to formula (2.10). 

Let's take some smooth orientable hypersurface a in M. By n = n(p) we de- 
note smooth field of unit normal vectors on a. Using Riemannian metric g in the 
manifold M, we can transform it to unit covector field: 

n 

ni = ^gijn^ . (2.11) 

We denote it by the same symbol n. Using covector n = n(p), we set up the 
following boundary- value problem for Hamilton- Jacobi equation (2.4): 

S\ =so = const, VSI = u{p) ■ n{p). (2.12) 

Here v = iy{p) is a scalar factor. It is not arbitrary. In order to make (2.12) 
compatible with the equation (2.4) it should satisfy the following condition: 

H{p,u(p)-n(p)) = 0. (2.13) 

Aside from boundary-value problem (2.12) for the equation (2.4), we set up the 
following initial-value problem for Hamilton equations (2.6): 

x'\ =x\p), pi\ =u{p)-ni{p). (2.14) 

lt=0 l*=0 

Initial data of Cauchy problem (2.14) are parametrized by points of hypersurface 
a. Function v = v{p) in (2.14) is the same scalar factor as in (2.12). It is restricted 
by condition (2.13) above. 

Solving Cauchy problem (1.18) for the equations (2.6), we get a family of para- 
metric curves in T*M. Their projections to M form a family of parametric curves 

l^l{t,p) (2.15) 

starting at the points of hypersurface a. Boundary data (2.12) and initial data 
(2.14) are called regular if function v = u{p) does not vanish on a, i. e. if 

u{p) ^ 0, (2.16) 

and if integral curves (2.15) are transversal to hypersurface a. Now, using Hamilton 
function H, let's define the following extended scalar field in M: 

n 

n = Y^PiV'H. (2.17) 

i=l 

In terms of transversality condition 7 ^ ct is expressed by inequality 



O(p,z/(p)-n(p))5^0. 



(2.18) 
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Suppose that both regularity conditions (2.16) and (2.18) on hypersurface cr are 
fulfilled. In this case integral curves (2.15) fill some neighborhood of hypersurface 
(T and we can formulate the following theorem. 

Theorem 2.1. If boundary data (2.12) satisfying condition (2.13) are regular, i. e. 
if both inequalities (2.16) and (2.18) on a are fulfilled, then boundary-value problem 
(2.12) for Hamilton- Jacobi equation (2.4) has unique solution S in some neighbor- 
hood of hypersurface a. This solution is given by integral 

t 

/n „ n 

dx' = So + j^Pix' dt. (2.19) 
^ i=i i i=i 

Theorem 2.1 is standard result in the theory of first order partial differential equa- 
tions. Its proof can be found in paper [1]. 

3. Wave front dynamics. 

Theorem 2.1 gives a method for solving Hamilton- Jacobi equation and finding 
scalar field S = S{p). Remember that S in (1.5) is a phase of propagating wave. 
By definition wave front is a set of points with constant phase: 

cr(s) = {pe M: ^(p) =s = const}. (3.1) 

Hypersurface a in (2.12) is interpreted as initial wave front, since cr = (t(.so). In 
regular case for s sufficiently close to sq wave fronts (j{s) in (3.1) are smooth hy- 
persurfaces filling some neighborhood of initial hypersurface a. Curves (2.15) are 
interpreted as rays in the limit of geometric optics (A oo), while Hamilton equa- 
tions (2.6) determine the motion of points, starting at the time instant t — Q and 

moving along these rays. If we fix their 
positions at some nonzero instant of time 
t, for sufficiently small value of t we would 
obtain some smooth hypersurface at close 
to initial hypersurface a. Hypersurfaces 
at form another family filling some neigh- 
borhood of initial wave front a. Should 
they coincide with a{s) ? In general, the 
answer is negative: 

at ^ a{s), (3.2) 

sec Fig. 3.1. Hypersurfaces at arc drawn 
in solid lines, while wave fronts a{s) are 
drawn in dashed lines. Inequality (3.2) 
means that Hamilton equations (2.6) do not describe real dynamics of wave fronts 
exactly. In order to describe this dynamics in terms of ODE's we should change 
parameter t in (2.6) for another parameter s = S{t). Formula (2.19) for scalar field 
S yields the relation of these two parameters: 

^ = J^p^i^ = Y,p^Y^H = n. (3.3) 



8 



R. A. SHARIPOV 



Changing t for s according to the formula (3.3), we keep symbol t for denoting this 
new parameter. Then from (2.6) we derive modified Hamilton equations 

where Q is given by formula (2.17). Differential equations (3.4) are called differ- 
ential equations of wave front dynamics. 

4. Legendre transformation. 

Hamilton equations (2.6) and modified Hamilton equations (3.4) are the equa- 
tions in cotangent bundle T*M . Now we shall transform them into the equations in 
tangent bundle TM. Note that for each point q = {p, p) of cotangent bundle T*M 
we can consider vector v e Tp{M) with components 

= W'H. (4.1) 

Joining together p and v, we get a point q = (p, v) of tangent bundle TM. This 
determines a map A""'^ : T*M TM, which is known as inverse Legendre trans- 
formation. We shall assume it to be invertible and denote by A direct Legendre 
transformation: A : TM — > T*M. Suppose that X is some extended tensor field. 
According to the definition 2.1, it is tensor- valued function in T*M. Then com- 
position X o A is tensor-valued function in tangent bundle TM. So we get another 
definition of extended tensor field. 

Definition 4.1. Extended tensor field X of type (r, s) in M is a tensor- valued 

function that maps each point q G G C TM to a tensor of the space T^{p,M), 
where p = 7r(g). Subset G of TM is a domain of extended tensor field X. If 
G = TM, then X is called global extended tensor field in M. 

In the case of arbitrary smooth manifold M definitions 2.1 and 4.1 lead to differ- 
ent theories. But for Ricmannian manifold M tangent bundle TM and cotangent 
bundle T*M are bound with each other by duality maps 

g : TM T*M, g-^ : T*M TM. (4.2) 

In local coordinates duality maps (4.2) arc represented as index lowering and index 
raising procedures in arguments of X. Due to duality maps (4.2) two objects 
introduced by definitions 2.1 and 4.1 are the same in essential. We call them 
covariant and contravariant representations of extended tensor field X. If we take 
contravariant representation of X, then formulas (2.9) and (2.10) are rewritten as 



dpi 



^.Xl:t - ^37^' (4-3) 



o=l 6=1 

r n s n 

j_ \ ' \ ' pifc y-il--- o-k-.-ir _ \ ' \ ' pfefc vil 

"^Z^ Z^^qa.k^3i Je 2^^ qjkJi-bk- 

fc=lafe = l fc=l 6^5=1 



(4.4) 
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Direct and inverse Legendre transformations are other two maps 

A : TM ^ T*M, A"^ : T*M TM (4.5) 

binding TM and T*M. Presence of nonlinear maps (4.5) increases the number 
of representations of extended tensor field X. If Y = Xo A~^, then we say that 
Y is p-representation or momentum representation for X, while X is called 
v-representation or velocity representation for Y. Total set of representations 
for extended tensor field is given in the following table: 



covariant p-representation 


contravariant p-representation 


covariant v-representation 


contravariant v-representation 



The following extended scalar field is called Lagrange function: 

n 

l = Y^Pi-V'H - H (4.6) 

Covariant p representation is natural for Hamilton function H, while Lagrange 

function is naturally used in contravariant v-rcpresentation: L = I " \. If Lagrange 
function (4.6) is already transformed to contravariant v-representation, then direct 
Legendre transformation A in (4.5) can be given by formula similar to (4.1): 

P^ = V^L (4.7) 

(see [3] for more details). Note that formulas (2.9), (2.10), (4.3), and (4.4) for 

V and V are valid either in p-representation and in v-representation. However, in 
v-representation we should replace pi by Vi and by v'^ in these formulas. Therefore 

V is called velocity gradient in v-representation. 

Duality maps (4.2) commutate with covariant differentiations V and V. This is 
expressed by the following equalities: 

V(Xog) = (VX)og, V(Xog) = (VX)og. (4.8) 

Legendre transformation A does not commutate with differentiations V and V. Un- 
like (4.8), we have the following equalities, where C is the operation of contraction: 

V(YoA) = C(VYoA(» WL), (4.9) 
V(YoA) = VYoA-t-C(VYoA(g) VVL). (4.10) 

In local coordinates the equalities (4.9) and (4.10) are written as follows: 

n 

%^t't = E ^^'^^ • ' (4-11) 

k=l 

n 

v^xi^::X = + ^ v^v^l ■ v'^y]^:::;;. (4.12) 

fe=i 
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Similarly, one can write inverse relationships for (4.9) and (4.10): 

V(XoA-i) =C(VXoA-i(8)VVif), (4.13) 
V(X o A"^) = VX o A"^ + C(VX o A"^ VVH). (4.14) 

In local coordinates these relationships (4.13) and (4.14) are expressed as 

n 

^kyii-ir- = V''V«if • V„X^-V, (4.15) 
q=l 

n 
9=1 

In (4.11), (4.12), (4.15), and (4.16) wc assume that X = YoA and Y = XoA'^. 

Now let's consider Lagrange function L. Its p-representation is given by formula 
(4.6). Let's calculate VfeZ directly by means of formula (4.6): 

n 

Vfe/ = Yj'i- "^k^'H - VkH. (4.17) 
Then let's apply formula (4.16) and let's calculate this derivative VfeZ again: 

n 
9=1 

Comparing (4.17) and (4.18), and taking into account (4.7), we get: 

VkL = -VkH. (4.19) 

In coordinate free form this equality (4.19) looks like VL = — ViJoA. 

5. Lagrangian representation for the 
equations of wave front dynamics. 

Now we are ready to transform differential equations (3.4) to v-representation. 
Using formula (4.1), we write fist equation (3.4) as 

Extended scalar field O in p-representation is determined by formula (2.17). Due 
to (4.1) and (4.7) its v-representation is determined by formula 

n 

n = Y^v' -S/iL. (5.2) 

i=l 

In order to transform second equation (3.4) we use formulas (4.7) and (4.19): 

Vt(V,L) = (5.3) 
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Written together, (5.1) and (5.3) form Lagrangian representation for the equations 
of wave front dynamics (3.4). Here are these equations: 

= V.(V.L) = ^. (5.4) 

Though written in terms of Lagrange function, the equations (5.4) are not Euler- 
Lagrange equations due to the denominator 17 in them. 

6. Dynamical systems admitting normal shift. 

Note that vector v with components v^, ... , i)" is not an actual velocity vector 
for the dynamics described by differential equations (5.4). Let's denote actual 
velocity vector by u. Then from first equation (5.4) we derive 

= -. (6.1) 

Denominator in (6.1) is not constant, it depends on a;-'^, ... , .t" and on compo- 
nents of velocity vector v as well. Therefore (6.1) defines nonlinear map 

H-.TM^TM (6.2) 

similar to Lcgcndrc map A determined by (4.7). If this map (6.2) invertible, then 
one can transform (5.4) to Newtonian form: 

x'=u\ Vtu' = F'(a;\... ,a;",u\... ,u"). (6.3) 

Components of force vector F in (6.3) are expressed through Lagrange function 

L and its derivatives. Now we shall not derive general formula for F''. But we 
consider very important special case, when L is fiberwise spherically symmetric 
extended scalar field in Riemannian manifold M. 

Definition 6.1. Extended tensor field X is called fiberwise spherically symmetric 
if its components depend only on modulus of velocity vector: 

^J:;;^: =^";;}:(a;\ ••• ,ar",t;), where v=\v\. 

This means that ^jlZjl is spherically symmetric function within each fiber of 
tangent bundle TM for each fixed point p G M. 

Fiberwise spherically symmetric tensor fields were considered in Chapter VII 
of thesis [4], see also recent papers [3] and [5]. Now suppose that L is fiberwise 
spherically symmetric extended scalar field: L = L{x^, ... , a;", w). Then from (5.2) 
we derive that O is also fiberwise spherically symmetric scalar field: 

n = |v| • L'. (6.4) 



Here L' is partial derivative of L with respect of its last argument v = |v|. Substi- 
tuting (6.4) into the equality (6.1), we derive the following formulas: 
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By A^' in (6.5) we denote components of unit vector N directed along vector v. 
Using its components, we define two operators of orthogonal projection Q and P: 

G-. = ^. n = 4-^. (6.6) 

First of them is a projector to the direction of vector v, second is a projector to 
hyperplane perpendicular to v. Projection operators Q and P with components 
(6.6) are complementary to each other, this means that 

Q + P = l and QoP = PoQ = o. 

Now let's consider second equation (5.4). For Vt(V,i) in it we have 

n n 

Vt(ViL) = J2 ^kViL ■x'' + Y, ^kViL ■ VtvK (6.7) 

fc=i fc=i 

By direct calculations for various derivatives in (6.7) we obtain 

WiL = L'-^=L'- Ni, (6.8) 

|v| 

VfcV,L = VfcL' • iV„ (6.9) 

VkViL = L" -Qik + ^-Pik. (6.10) 
|v| 

Substituting (6.9) and (6.10) into (6.7), we get the following equality: 

VtiViL) = ^ ^ . g^fc + ^ ( L" . Q,, + . P,, . V^^;^ (6.11) 

k=l k=l ^ ^ 

Next step is to relate Vtv'^ in (6.11) with Vtu''. For this purpose let's apply Vt to 
the first equality in (6.5). As a result of direct calculations we get: 

V.«' = t ( ■ - |p - Qi) ■ V.»- - g ^ - (6.12) 

In order to express VtV^ through VtU*^ we should invert matrix D with components 

1 T" 

P)k pfe _ f\k 

' \v\-L' ' {L'f'^'- 

This matrix is invertible if and only if L" ^ 0. The condition L' 7^ is fulfilled 
since = |v| • L' 7^ 0. This follows from transversality condition (2.18). Inverting 
matrix D and applying it to (6.12), we get the following expression for VtV*: 

V,.'^ = (|v| . L' . Pi= - . Ql) ■ V,u^ - E ^ • Q''- (6.13) 

s=l ^ ' s=l 
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Now let's substitute (6.13) back to right hand side of (6.11). As a result we get 

n 

VtiViL) = {L'f Y,{Pik - Qik) ■ Vtu". (6.14) 

fe=i 

Then let's substitute (6.14) into second equation (5.4) and let's use formula (6.4) 
for denominator Cl in it. This leads to the equation for VjU*': 

" X7 r 

[L'f Y,{Pik - Qik) ■ Vtu" = -f--. (6.15) 
Equation (6.15) can be explicitly solved. For Vtu'' from (6.15) we derive 

" \7 r 

Vt«' = E ^T(I7)3 • ia'' -^N'N')- (6.16) 

Comparing (6.16) with (6.3) we see that right hand side of (6.16) is the expression 
for F'^. Formula for covariant components is more elegant: 

" V r 

The only problem now is that VjZ/, L', and |v| in right hand side of (6.17) are 
given in v-representation, while components of force vector in (6.3) should de- 
pend on components of actual velocity u. Our last effort is to transform (6.17) to 
u-representation, using map (6.2). 

Let's denote by h Hamilton function H in v-representation, i.e. h = H o X. Then 
for h we have the following formula similar to formula (4.6) for I: 

n 

h = ^v'-ViL-L=\v\-L'-L. (6.18) 

i=l 

Differentiating (6.18), we obtain 

Vih = v ViL' - ViL, h' = v- L". (6.19) 

Let W be u-rcprcscntation for h. Denote by u modulus of actual velocity vector. 
From (6.5) and (6.18) we conclude that W is & fiberwise spherically symmetric 
scalar field in u-representation: W = W{x^, . . . ,a;",u), where u = |u|. Now let's 
calculate derivatives of W. For W and ViW we get: 

du dv du \9v J ' 

(6.20) 

dW . , du ( du\~^ 
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Now let's remember that u and v are related by second equality in (6.5). We write 
it as « = e/L' , where e = ±1 depending on sign of L' . Then we get 



du e-L" du £-ViL' 



(6.21) 



dv {L'Y' dx' {L'Y 

Combining (6.19), (6.20), and (6.21), we can transform (6.22) as follows: 

W' = -e- |v| • [L'f, WiW = -ViL. (6.22) 

If we write second equality (6.5) as |u| = e/L', then from (6.22) we derive: 



W |v| • {L'Y 



(6.23) 



Comparing the equality (6.23) with formula (6.17) for components of force vector 
F, we can rewrite this formula in the following form: 

Fk = -|u| .j2^-{2N'N,-Si). (6.24) 

1=1 

Now formula (6.24) is completely compatible with (6.3). Its right hand side is given 
in u-representation as it is required in (6.3). 

Formula determines force field of Newtonian dynamical system (6.3) describing 
the dynamics of wave front for linear wave equation (1.4). But most important fact 
for us is that formula (6.24) gives the link to the theory of Newtonian dynamical 
systems admitting normal shift of hypersurfaces. This theory was developed in 
series of papers [6-21], on the base of which theses [4] and [22] were prepared. 
According to results of Chapter VII of thesis [4] (see also paper [5]), force field 
of any Newtonian dynamical system admitting normal shift of hypersurfaces in 
Riemannian manifold of the dimension n 3 is given by explicit formula: 

^. = MEl^-|u|.t^.(2*<*-«. (6.25, 

Here W = W{x^, ... , a:", |u|) is arbitrary fiberwise spherically symmetric extended 
scalar field with W ^ 0, and h = h{w) is arbitrary function of one variable. 

7. Conclusions. 

Comparing formulas (6.24) and (6.25), we see that they do coincide for the case 
h = 0. Term with h = h{W) in (6.25) is responsible for energy dissipation and 
energy pumping phenomena. In the absence of these phenomena, i. e. for h = 0,we 
can state the following main results of present paper: 

1) conservative dynamical systems admitting normal shift of hypersurfaces in Rie- 
mannian manifold of the dimension n ^ 3 coincide with those describing wave 
front dynamics in quasiclassical limit for wave operators with fiberwise spheri- 
cally symmetric symbol H; 

2) these systems are not Hamiltonian, but they are very close to Hamiltonian 
systems (see equations (3.4) above); 

3) they possess first integral, which can be interpreted as energy (or as Hamilton 
function in p-representation) . 
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Third result that is a first integral of Newtonian dynamical system with force 
field (6.24) was already known (see paper [5] or Chapter VII of thesis [4]). But its 
interpretation was not so clear as it is now. 

Note that the equations of wave front dynamics do not exhaust the whole class of 
Newtonian dynamical systems admitting normal shift of hypersurfaces in M . The 
problem of proper interpretation for non-conservative term 



KW) ■ Nk 

w 

in formula (6.25) in the sense of theory from [1] and (23) is still open. This problem 
should be studied in separate paper. Methods developed in book [24] might appear 
useful for this purpose. 

Remark. Technique based on usage of extended tensor fields is not new (see for 
instance book [25]). And, as noted by Prof. Mircea Crasmareanu, our terminology 
above is not standard. However, we do not change this terminology in order to 
keep integrity of present paper and our previous papers (see [3-22] and some other 
papers in Electronic Archive at LANL^). 
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